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Diffusion has been widely used to describe a random walk of particles or waves, and it re- 
quires only one parameter - the diffusion constant. For waves, however, diffusion is an ap- 
proximation that disregards the possibility of interference 1 . Anderson localization 2 , which 
manifests itself through a vanishing diffusion coefficient in an infinite systemPE originates 
from constructive interference of waves traveling in loop trajectories - pairs of time-reversed 
paths returning to the same point 5 6 . In an open system of finite size, the return probabil- 
ity through such paths is reduced, particularly near the boundary where waves may escape. 
Based on this argument, the self-consistent theory of localization and the supersymmetric 
field theory predict 7 9 that the diffusion coefficient varies spatially inside the system. A di- 
rect experimental observation of this effect is a challenge because it requires monitoring wave 
transport inside the system. Here, we fabricate two-dimensional photonic random media and 
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probe position-dependent diffusion inside the sample from the third dimension. By varying 
the geometry of the system or the dissipation which also limits the size of loop trajectories, 
we are able to control the renormalization of the diffusion coefficient. This work shows the 
possibility of manipulating diffusion via the interplay of localization and dissipation. 

As first shown by Einstein in his theory of Brownian motion, the diffusion equation describes 
the evolution of the density of particles each undergoing a random walkP^. The concept of diffu- 
sion has since been used to describe transport phenomena in both physics and other sciences. The 
power of this approach is that it requires knowledge of a single parameter, the diffusion constant, 
regardless of the underlying microscopical mechanisms of transport. If the spatial gradient of 
particle density is not too large, the particle flux is linearly proportional to the gradient, and the 
coefficient is the diffusion constant. Diffusion is also applicable to waves-^l but it ignores inter- 
ference effects. When inelastic scattering is negligible, most of the elastically scattered waves have 
uncorrelated phases and their interference is averaged out. Nevertheless, a wave may return to a 
position it has previously visited after a random walk, and there is always the time-reversed path 
which yields identical phase delay. Constructive interference of the waves from the reversed loops 
increases wave (energy) density at the original position and decreases the flux, giving the so-called 
weak localization effecP. This is the basic mechanism for the suppression of wave diffusion, which 
eventually leads to Anderson localization^. 

In the self-consistent theory of localization, the diffusion coefficient D is renormalized, and 
the amount of renormalization is proportional to the return probability of waves via the looped 
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paths 3 4 . In an open system of finite size, the return probability is reduced because the longer loops 
may reach the boundary where waves escape. Thus the renormalization of D depends on the sys- 
tem size. Moreover, near the boundary the chance of escape is higher, so the renormalization of D 
is weaker. This means the value of D is no longer constant but varies spatially^. In the presence of 
dissipation the long loops are also cut, thus the renormalization of the diffusion coefficient depends 
on the amount of dissipation. This sets an effective system size beyond which the wave will not 
return 15 . Although self-consistent theory has successfully interpreted several experiments^^, its 
key prediction of position-dependent diffusion has not been observed directly because it is difficult 
to probe wave transport inside the system experimentally. 

Here we report direct experimental evidence of position-dependent diffusion by probing light 
transport inside a quasi-two-dimensional random system from the third dimension. The system 
size and shape are designed to make the return probability sufficiently high so that the diffusion 
coefficient is modified appreciably. We also use dissipation to control the effective system size, 
and tune the value of D via the interplay of localization and dissipation. This work demonstrates 
the possibility of utilizing the geometry of a random system or the dissipation to manipulate wave 
diffusion. 

We designed and fabricated two-dimensional (2D) disordered waveguide structures in a 
220 nm silicon layer on top of 3 /im buried oxide. The patterns were written by electron beam 
lithography and etched in an inductively coupled plasma reactive ion etcher. As shown in the scan- 
ning electron microscope images in Fig. [TJ the waveguide has sidewalls made of periodic arrays 
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of air holes. They possess a 2D photonic bandgap that covers the wavelength range of the probe 
light (A = 1500 nm - 1520 nm), thus providing optical confinement in the plane of the waveg- 
uide. Light enters the waveguide from an open end and is incident onto a 2D array of air holes 
inside the waveguide. The random pattern of air holes causes light to scatter while going through 
the waveguide. The transport mean free path t is determined by the size and density of air holes. 
Light localization will occur if the length of the random array L exceeds the localization length 
£ = (tt/2)N£, where iV = 2W/ (A/n e ) is the number of propagating modes in the waveguide, W 
is the waveguide width, A is the optical wavelength in vacuum, and n e is the effective index of 
refraction of the random medium. Since N scales linearly with W, £ can be easily tuned by vary- 
ing the waveguide width. Therefore, by changing the waveguide geometry (L, W), we can reach 
both the diffusion regime (£ < L < £) and localization regime (L > £j 20 * 21 [ Although there is no 
mobility edgepl in such a system, it is not essential for our goal of observing position-dependent 
diffusion. 

In order to apply the self-consistent theory of localization to the analysis of the experimen- 
tal data below, we first validate it with numerical simulations under conditions close to those in 
the experiment. We computed the position-dependent diffusion coefficient without making any 
assumption about the nature or strength of wave interference (see the Methods section). Figure |2] 
plots the calculated D(z)/D (the lower dashed line) for L/£ = 3.0, where D is the diffusion 
coefficient without renormalization.. The z axis is parallel to the waveguide, and the random ar- 
ray extends from z = to z = L. The renormalized D(z) drops to 0.17D in the middle of the 
random waveguide (z = L/2). Using the self-consistent theory of localization (see the Supple- 

4 



mentary Information) we calculate D(z)/D (the lower solid line in Fig. |2) and it is in excellent 
agreement with the ab-initio simulation without any fitting parameters. Previous studies show that 
further into the localization regime where resonant tunneling dominates wave transport, the self- 
consistent theory of localization underestimates the energy density inside the random system that 
is strongly affected by the presence of necklace states^. In our experiment we avoid such a regime 
and stay where the self-consistent theory of localization holds^. 

Another factor we shall consider is the dissipation of light in the random waveguide. Within 
the wavelength range of the probe, light absorption by silicon or silica is negligible. However, light 
is scattered out of the waveguide plane by the random array of air holes. Such scattering allows us 
to monitor the intensity distribution inside the system from the vertical direction. The question is 
whether the out-of-plane scattering can be treated as incoherent dissipation. In a periodic array of 
scatterers, long-range correlation of light fields makes the waves scattered from different locations 
phase coherent and their interference in the far field zone determines the out-of-plane leakage. In 
contrast, in a random array of scatterers, the fields are correlate d 24 * 25 ^ only within a distance of 
the order one transport mean free path £, and waves from different coherent regions of size I x I 
have uncorrected phases. Since there are a large number of such coherence regions t x i in our 
waveguides W x L, the overall leakage may be considered incoherent and treated effectively as 
absorption. 

To illustrate the effect of dissipation on position-dependent diffusion, we perform numerical 
simulations. The diffusive absorption length in the random system is £ a0 = \fDtfF a , where r a is 
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the ballistic absorption time. When £ a0 becomes smaller than the localization length £, the effect 
of dissipation is significant. Figure |2] plots the calculated D(z)/D in the random waveguide 
with £ a o/£ = 0.45 (the upper dashed line) in comparison to that with £ a0 = oo (no absorption). 
The suppression of diffusion is weakened by the absorption, and a plateau for the renormalized 
diffusion coefficient is developed inside the disordered system. This result can be understood as 
follows. Dissipation suppresses feedback from long propagation paths, limiting the effective size 
of the system 21 to the order of diffusive absorption length for any position that is more than one 
£ a0 away from the open boundary (£ a0 < z < L — £ a o) 23 - Thus the renormalized D reaches a 
constant value equal to that of an open system of dimension ~ 2£ a = 2\/Dr a . In the remaining 
regions that are within one £ a o to the boundary (z < £ a o and L — z < £ a0 ), the diffusion coefficient 
is still position dependent due to leakage through the boundary and D increases toward the value 
of D without renormalization. We note that the extent of these regions £ a0 is much greater than 
the transport mean free path £. The latter determines the boundary region where the diffusion 
approximation is not accurate even without wave interference^. Figure |2]also shows the prediction 
of self-consistent theory of localization in the presence of dissipation (the upper solid line), and it 
agrees well with the numerical result. 

Experimentally there are three advantages to using the planar waveguide geometry. First, it 
allows a precise fabrication of the desired system using lithography so that the parameters such as 
the transport mean free path can be accurately controlled. Second, the localization length £ oc W 
can be varied by changing the waveguide width, while the diffusive absorption length £ a0 remains 
fixed. This allows us to separate the effects of localization and dissipation by testing waveguides 
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Figure 1 : Experimental realization of disordered waveguides. Top-view a and tilt- view b scan- 
ning electron microscopy (SEM) images of an optical waveguide fabricated in a silicon membrane 
on top of silica. The two sidewalls of the waveguide consist of triangular lattices of air holes (lat- 
tice constant 440 nm, hole radius 154 nm). They possess a 2D photonic bandgap and behave like 
reflecting walls for light incident from all angles in the waveguide. The probe light is coupled from 
a silicon ridge waveguide to an empty photonic crystal waveguide, then impinge onto a random 
array of air holes (hole diameter 100 nm, and areal density 6 %) inside the waveguide. 
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of different width. Finally, unlike 2D random systems 27 , the additional confinement of light by the 
waveguide sidewalls makes £ scale linearly with I. Even if scattering is relatively weak (k£ 3> 1, 
where k is the wavenumber), the waveguide length L can easily exceed £ so that the localization 
effect is strong enough to modify diffusion. Instead of designing the disorder to maximize scat- 
tering (minimizing k£), we deliberately lower the density of air holes to mitigate the out-of-plane 
scattering loss and maximize the ratio £ a o/£. 

In the optical measurement, output from a wavelength-tunable continuous-wave laser (HP 
8168F) was coupled to the waveguide through a single-mode polarization-maintaining lensed fiber 
as shown schematically in Fig. [3^. To ensure efficient confinement in disordered waveguide, the 
transverse-electric (TE) polarization (electric field in the plane of the waveguide) of the incident 
light was chosen. A near-field optical image of the spatial distribution of light intensity across 
the structure surface was taken by collecting light scattered out of plane using a 50X objective 
lens (numerical aperture 0.42) and recorded by an InGaAs camera (Xenics Xeva 1.7-320). An 
example of the near-field image is shown in Fig. |3b, from which the intensity distribution I(y, z) 
is extracted. The y axis is in the waveguide plane and perpendicular to the direction of incident 
beam. I(y, z) was integrated over the cross section of the waveguide (along the y axis) to obtain the 
evolution of intensity I(z) in the incident direction (parallel to the z axis). For each configuration 
(width W, length L, transport mean free path £) of the disordered waveguides, I(z) was averaged 
over two random realizations of air holes and fifty input wavelengths equally spaced between 
1500 nm and 1520 nm. The wavelength spacing was chosen to produce independent intensity 
distributions. 
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Figure 2: Numerical simulations of position-dependent diffusion. Numerically calculated 
diffusion coefficient D(z) (dashed curves) at position z (0 < z < L) in a waveguide filled with 
randomly positioned scatterers. The waveguide has length L = 3£ (£ is the localization length) and 
supports N = 10 modes. The solid curves represent the prediction of the self -consistent theory of 
localization. D denotes the diffusion coefficient that ignores interference effects. In the absence 
of absorption (diffusive absorption length £ a0 = oo), D(z) drops to a minimum of 0.17D in the 
middle of the waveguide. With the addition of absorption (£ a o/£ = 0.45), D(z) exhibits a plateau 
for £ a0 < z < L — £ a0 , and its value D p is determined by the ratio £ a o/£- 
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Figure 3: Optical measurement of light inside random medium, a, A schematic depiction of 
light coupling from a lensed fiber to a waveguide in the silicon membrane, b, A near-field optical 
image showing the spatial distribution of light intensity in a random waveguide with a configuration 
similar to that shown in Fig. [Q The wavelength of the probe light is 1510nm. A50x objective lens 
collects the light scattered by the air holes out of the waveguide plane and images onto a camera. 
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Figured shows the measured I(z) inside the ensembles of random waveguides of W varying 
from 60 fim to 5 /im (blue solid lines). All other parameters are kept the same, L is fixed at 80 /mi, 
the diameter of air holes is 100 nm, and average (center-to-center) distance of adjacent holes is 
390 nm. £ and £ a0 are obtained by fitting the least localized sample, W=60 /mi (longest £), with 
the self-consistent theory of localization (red dashed line). We find that £ = 2.2 ± 0.1 fim and 
£ a o = 30 ± 0.5 /im. With the parameters found from the 1^=60 /im sample, the self-consistent 
theory of localization successfully predicts the decay for I(z) in all other samples with 1^=40 /mi, 
20 /im, 10 /im, 5 /im (red solid lines). We stress that the excellent agreement with the experimental 
data are obtained without any free parameter except the vertical intensity scale. The position- 
dependent diffusion coefficients D(z) corresponding to the red curves in Fig. @^ are shown in 
Fig.BJ). We can clearly see that the diffusion coefficient is reduced inside the sample, and its value 
varies along z. Farther away from the open boundary, D has a smaller value. In the narrower 
waveguides, the reduction of D is larger due to stronger localization effect. In the most localized 
sample of W = 5/mi, D is reduced to 0.65-D at 2 = L/2. In an attempt to further reduce D, we 
double the length of random system L to 160 /im. As shown in Fig.|4fc,d for W=5 /im, the minimal 
D no longer decreases, instead it saturates in the middle of the random waveguide. This behavior 
is attributed to dissipation which suppresses localization. As the system length L becomes much 
larger than the diffusive absorption length £ a0 , D(z) saturates to a constant value D p inside the 
disordered waveguide, similar to the simulation result shown in Fig. [2] 

Finally we exploit the interplay between dissipation and localization to tune the saturated 
value of the diffusion coefficient inside the random system. To this end, we increase the density of 
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Figure 4: Experimental measurement of light transport inside disordered waveguides, a, Ex- 
perimentally measured light intensity I(z) inside random waveguides of different width W and 
constant waveguide length L=80 /im (blue solid lines). The curves are vertically shifted for a 
clear view. £=2.2 /mi and £ a o=30 /tm are found by fitting the 1^=60 /im sample with the self- 
consistent theory of localization (red dashed line). With these parameters, the self-consistent the- 
ory of localization predicts I(z) for other samples of 14^=40 /im, 20 /tm, 10 /im, 5 /tm (red solid 
curves), and the prediction agrees well with the experimental data without any free parameter ex- 
cept the vertical intensity scale, b, Position-dependent diffusion coefficients for the five samples 
in a. c, Experimentally measured I(z) of two waveguides with the same width W=5 /im but dif- 
ferent length, L=80 /im, 160 /tm (blue solid curves). Red solid curves represent the prediction of 
the self-consistent theory of localization using the same values of £ and £ a0 as in a. d, Diffusion 
coefficients D(z) for the two samples in c, showing the saturation of D inside the longer sample 
L=160 /tm. 



scatterers to reach the deep saturation region £ a0 <C L. In the second set of samples, the diameter of 
air holes is 150 nm, and the average distance between adjacent holes is 370 nm. Waveguide length 
L is set at 80 /im and W varies from 5 /im to 60 /im. Experimental data of measured intensity 
I(z) inside the random waveguides are presented in Fig. |5Jl. By fitting the W=60 /im sample 
with the self-consistent theory of localization (red dashed line), we obtain i = 1.0 ± 0.1 /im 
and £ a0 = 13.0 ± 0.3 /mi. With these values, the self-consistent theory of localization gives the 
spatial profiles of I(z) for 1^=40 /zm, 20 /mi, 10 /im, 5 /zm (red solid lines), which agree well 
with the experimental data with no fit. The corresponding values of D(z) are plotted in Fig. |5b. 
Due to stronger scattering (smaller t) and larger out-of-plane loss (shorter £ a0 ), D(z) for all five 
samples displays a well-developed plateau inside the sample. The saturated value D p decreases 
with W - the narrower waveguide has smaller D p . Hence, without changing the disorder or altering 
the dissipation rate, we can control the diffusion inside a random system by merely varying its 
geometry (W in this case). 

Renormalization of the diffusion coefficient D has long been considered as a theoretical 
approach to tackle the problem of localization. Here we present direct experimental evidence of 
suppressed diffusion of light inside the random systems. This is an intrinsic wave phenomenon, 
hence our conclusions also apply to acoustic wave, microwave and even the de Broglie wave of 
electrons. The diffusion coefficient D is altered by wave interference and such modification de- 
pends on the position inside an open random system. By varying the size and shape of the random 
system, we are able to control the strength of wave interference and consequently the degree of 
renormalization of D. We also show that the presence of dissipation prevents D from approach- 
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Figure 5: Tuning the diffusion coefficient via the interplay of localization and dissipation. 



a, Experimentally measured light intensity I(z) inside random waveguides in the deep saturation 

regime £ a0 <C L (blue solid lines). The curves are vertically offset for a clear view. The length and 

width of the waveguides are given in the graph. (. = 1.0 fim and £ a o = 13 fim are found by fitting 

the 1^=60 fim sample with the self-consistent theory of localization (red dashed line). This is then 

used to predict I(z) for other samples 1^=40 /zm, 20 /mi, 10 /jm, 5 fxm (red solid curves). The 

prediction of the self-consistent theory of localization is in good agreement with the experimental 

data with no fitting parameters except the vertical intensity scale, b, Diffusion coefficients D(z) 

for all samples in a are saturated in the region £ a0 < z < L — £ a0 . The saturated values are smaller 
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ing zero and sets a limit for the minimal D that can be reached by the localization effect. Such 
effect of dissipation is expected to be similar to that of dephasing in the electronic systems 28 . The 
interplay between localization and dissipation enables us to tune the value of D inside random 
systems. The results presented in this work are obtained by directly probing light transport inside 
the random media, which allows us to extract the diffusion coefficient anywhere inside the system. 
Such experiments open the path to measure other transport properties inside random systems, e.g., 
intensity correlations and fluctuations. 

Methods 

Numerical simulations: The numerical data shown in Fig. |2] were obtained from the ab-initio sim- 
ulation of a monochromatic scalar wave propagating in a 2D waveguide filled with random scat- 
terers. With continuous- wave excitation from one end of the waveguide (z = 0), we computed 
the energy density W(z) and flux along the z axis J z (z), and averaged them over the cross sec- 
tion of the waveguide. The position-dependent diffusion coefficient was found from Fick's law as 
D(z) = — ( J z (z)) / [d(W(z))/dz], where angular brackets denote the ensemble average. Using 
a supercomputer, we numerically simulated an ensemble of 10 6 waveguides with different disor- 
der configurations. To obtain the value of the diffusion coefficient D without renormalization, 
we used the single parameter scaling and the expression for direction-resolved flux, as detailed in 
Supplementary Information. In the absorbing samples we used the continuity equation to calculate 
the absorption time r a . 

Self-consistent theory of localization: The application of self-consistent theory of localization with 
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a position-dependent diffusion coefficient to disordered waveguides was described in Ref. It 
involves the diffusion equation which defines the return probability and a self-consistency equation 
that relates diffusion coefficient D(z) to the return probability. We solved these two equations by 
iteration until we found D(z) which satisfied both equations, see the Supplementary Information 
for more details. 

Design of photonic crystal walls for 2D waveguides: The triangular lattice of air holes that form the 
sidewalls of the random waveguide were designed to have a 2D photonic bandgap for TE polarized 
light in the wavelength range of 1450 nm - 1550 nm. The photonic band structure was calculated 
with the plane wave expansion method^. 

References 

1. van Rossum, M. C. & Nieuwenhuizen, T. M. Multiple scattering of classical waves: mi- 
croscopy, mesoscopy, and diffusion. Rev. Mod. Phys. 71, 313-371 (1999). 

2. Lagendijk, A., van Tiggelen, B. & Wiersma, D. S. Fifty years of anderson localization. Phys. 
Today 62, 24-29 (2009). 

3. Vollhardt, D. & Wolfle, R Diagrammatic, self-consistent treatment of the anderson localization 
problem in d < 2 dimensions. Phys. Rev. B 22, 4666-4679 (1980). 

4. Kroha, J., Soukoulis, C. M. & Wolfle, R Localization of classical waves in a random medium: 
A self-consistent theory. Phys. Rev. B 47, 1 1093-1 1096 (1993). 



16 



5. Gor'kov, L., Larkin, A. & Khmel'nitskii, D. Particle conductivity in a two-dimensional random 
potential. JETP Lett. 30, 228-232 (1979). 

6. Hikami, S. Anderson localization in a nonlinear-cr-model representation. Phys. Rev. B 24, 
2671-2679(1981). 

7. van Tiggelen, B. A., Lagendijk, A. & Wiersma, D. S. Reflection and transmission of waves 
near the localization threshold. Phys. Rev. Lett. 84, 4333-4336 (2000). 

8. Cherroret, N. & Skipetrov, S. E. Microscopic derivation of self-consistent equations of anderson 
localization in a disordered medium of finite size. Phys. Rev. E 77, 046608 (2008). 

9. Tian, C. Supersymmetric field theory of local light diffusion in semi-infinite media. Phys. Rev. B 
77, 064205 (2008). 

10. Einstein, A. liber die von der molekularkinetischen theorie der warme geforderte bewegung 
von in ruhenden fliissigkeiten suspendierten teilchen. Annalen der Physik 17, 549-560 (1905). 

11. Chandresekhar, S. Radiative transfer (Dover, New York, 1960). 

12. Ishimaru, A. Wave Propagation and Scattering in Random Media (Academic Press, 1978). 

13. van de Hulst, H. C. Multiple light scattering (Academic, New York, 1980). 

14. Anderson, P. W. Absence of diffusion in certain random lattices. Phys. Rev. 109, 1492-1505 
(1958). 

15. Payne, B., Yamilov, A. & Skipetrov, S. E. Anderson localization as position-dependent diffu- 
sion in disordered waveguides. Phys. Rev. B 82, 024205 (2010). 

17 



16. Storzer, M., Gross, P., Aegerter, C. & Maret, G. Observation of the critical regime near 
anderson localization of light. Phys. Rev. Lett. 96, 063904 (2006). 

17. Hu, H., Strybulevych, A., Page, J. H., Skipetrov, S. E. & van Tiggelen, B. A. Localization of 
ultrasound in a three-dimensional elastic network. Nat. Phys. 4, 945-948 (2008). 

18. Zhang, Z. Q., Chabanov, A. A., Cheung, S. K., Wong, C. H. & Genack, A. Z. Dynamics of 
localized waves: Pulsed microwave transmissions in quasi-one-dimensional media. Phys. Rev. B 
79, 063904 (2009). 

19. Sperling, T., Biihrer, W., Aegerter, C. M. & Maret, G. Direct determination of the transition 
to localization of light in three dimensions. Nat. Phot. 7, 48-52 (2013). 

20. Beenakker, C. W. Random-matrix theory of quantum transport. Rev. Mod. Phys. 69, 731-808 
(1997). 

21. Brouwer, P. W. Transmission through a many-channel random waveguide with absorption. 
Phys. Rev. B 57, 10526-10536 (1998). 

22. Tian, C., Cheung, S. & Zhang, Z. Local diffusion theory for localized waves in open media. 
Phys. Rev. Lett. 105, 263905 (2010). 

23. Yamilov, A. & Payne, B. Interplay of absorption and localization in absorbing disordered 
waveguides. arxiv:1303. 271 7 (2013). 

24. Shapiro, B. Large intensity fluctuations for wave propagation in random media. Phys. Rev. 
Lett. 57,2168-2171 (1986). 

18 



25. Feng, S., Kane, C, Lee, P. A. & Stone, A. D. Correlations and fluctuations of coherent wave 
transmission through disordered media. Phys. Rev. Lett. 61, 834-837 (1988). 

26. Akkermans, E. & Montambaux, G. Mesoscopic Physics of Electrons and Photons (Cambridge 
University Press, Cambridge, UK, 2007). 

27. Riboli, F. et al. Anderson localization of near-visible light in two dimensions. Optics Letters 
36, 127-129 (2011). 

28. Altshuler, B. L., Lee, P. A. & Webb, R. A. (eds.) Mesoscopic Phenomena in Solids (North 
Holland, Amsterdam, 1991). 

29. Johnson, S. G. & Joannopoulos, J. D. Block-iterative frequency-domain methods for 
maxwell's equations in a planewave basis. Opt. Express 8, 173-190 (2001). 

30. Mirlin, A. Statistics of energy levels and eigen-functions in disordered systems. Phys. Rep. 
326, 259-382 (2000). 

Acknowledgements We are indebted to Patrick Sebbah and Shivakiran Bhaktha for their insight in select- 
ing the experimental geometry. We acknowledge Seng Fatt Liew, Douglas Stone, Arthur Goetschy, Boris 
Shapiro, Chushun Tian and Sergey Skipetrov for useful discussions. We also thank Michael Rooks for 
suggestions regarding sample fabrication. This work was supported by National Science Foundation under 
grants Nos. DMR- 1205307, DMR-1205223 and ECCS-1 128542. Computational resources were provided 
under the Extreme Science and Engineering Discovery Environment (XSEDE) grant No. DMR-100030. 
Facilities use was supported by YINQE and NSF MRSEC DMR-1 119826. 



19 



Author contributions A.Y. and H.C. initiated the study and designed the experiments. R.S. and B.R. 
fabricated the samples. H.N. designed the waveguide walls. R.S. collected the experimental data. R.S. and 
A.Y. analyzed the data. B.R and A.Y. performed ab-initio numerical simulation. A.Y. and H.C. prepared the 
manuscript. 

Additional information Supplementary information is available in the online version of the paper. Reprints 
and permissions information is available online at www.nature.com/reprints. Correspondence and requests 
for materials should be addressed to A.Y. and H.C. 

Competing Interests The authors declare that they have no competing financial interests. 



20 



SUPPLEMENTARY INFORMATION 



Position-dependent diffusion of light in disordered waveguides 

Alexey G. Yamilov, Raktim Sarma, Brandon Redding, Ben Payne, Heeso Noh & Hui Cao 

1 Calculation of position-dependent diffusion coefficient D{z) 

In the ab-initio numerical simulation, we consider a monochromatic scalar wave _E(r)e~* w * prop- 
agating in a 2D volume-disordered waveguide of width W and length L ^> W. The wave field 
E(r) obeys the 2D Helmholtz equation: 

{V 2 + fc 2 [l + 5e(r)]}£(r) = 0. (1) 

Here k = to/c is the wavenumber and 6e(r) = (1 + ia)Se r (r), where 8e r (r) describes the random 
fluctuation of the dielectric constant, and a > denotes the strength of absorption. The system 
is excited from one open end (z = 0) of the waveguide (extending from z = to z = L) by 
illuminating each of the guided modes with a unit flux. The wave field E(r) throughout the random 
medium is computed with the transfer matrix method for a given realization of disorder 15 . From 
E(r) we calculate the energy density W(z) and the flux J z (z) along the z axis (parallel to the 
waveguide axis). These two quantities are averaged over the cross section of the waveguide at each 
z and give the diffusion coefficient: 

D(z) = -(J z (z))/[d(W(z))/dz], (2) 
1 



where the averages (. . .) are taken over a statistical ensemble of 10 6 disorder realizations. 

In order to compare our numerical results for D(z) with the self-consistent theory of local- 
ization, we need to have the value of the diffusion coefficient without renormalization due to the 
wave interference effects D = v£/2. To estimate the transport mean free path £ in our model we 
perform a set of simulations for different waveguide lengths L, exploring both the regime of diffu- 
sion L < £ and that of Anderson localization L > £. We computed numerically the conductance 
g as the sum of transmission coefficients from all incoming to all outgoing waveguide modes. The 
dependencies of the average (g) and variance var(g) on L are fitted by the analytical expressions 
obtained by Mirlin in Ref. 30 using the supersymmetry approach with i being the only fit parameter. 
To find the diffusive speed v we use the definition of diffusive flux in the forward (+z) direction 
J^ + \z) and the backward {—z) direction Jjf~\z) with respect to the propagation direction 1 

<jW(z)> = (v/-k)(W(z)) =F (D(z)/2)d(W(z)}/dz. (3) 

Combining the two components we find the diffusive speed 

v =2(<J<+>(z)) + <•/<->(*)>) /W*)>. W 
Dashed lines in Fig.[2]depict D(z) found in equation © normalized by D . 

In the dissipative random waveguides, the characteristic absorption time r a is determined 
numerically using the condition of flux continuity d (J z (z)) jdz = (l/r ) (W(z)). The desired 
diffusive absorption length £ a0 = \JTJtfF a can be obtained by the proper choice of a in equation (OQ). 
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2 Self-consistent theory of localization 



The self-consistent theory starts with the Green's function G(r, r') of equation (0Q) with 5e(r) = 
5e r (r)+ia. In a random waveguide, the disorder-averaged function C(r, r') = (4nWD /cL)(\G(r,r')\ 2 ) 
obeys self-consistent equations in a dimensionless formP^: 



' l\ 2 d Jf ^d' 



c(£0 = *(C-0, (5) 



1 or, A 

1 + — C(CC), (6) 



d(0 e 

where d(Q = D{Q/Dq and all position-dependent quantities are functions of the longitudinal 
coordinate ( = z/L. The quantity C((, (), which renormalizes the diffusion coefficient, is pro- 
portional to the return probability at £. Assuming first that d(() = 1, equations (15161) are solved by 
iteration with the boundary conditions: 

c(£ 0*^(0^6(6 O = o (7) 

at C = and £ = 1, The zo = (tt/4)£ is the so-called extrapolation length 1 . 

After the self-consistent solution of equations d5]|7]) has been found, we find the inten- 
sity distribution inside the sample by replacing the delta- function source in equation © with 
(L/£)exp [—(/(£/ L)]. This source term represents the exponential attenuation of the incident 
ballistic signal. 
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3 Experimental details 



The experimental setup for optical characterization is shown in Fig. SI (a). We used a single-mode 
polarization-maintaining fiber to deliver the probe light into a silicon ridge waveguide on a SOI 
substrate. The fiber was tapered at the end to focus the laser beam to a spot of diameter ~ 2.5 /zm 
at the edge of the wafer. The ridge waveguide had the same width as the random waveguide it 
was connected to, which varied from 5 micron to 60 micron [Fig. Sl(b)]. However, the height 
of the silicon waveguide was merely 220 nm, so some of the input light did not couple into the 
waveguide; instead it propagated above or below the waveguide. To avoid such stray light, the 
ridge waveguide was tilted by 30 degrees with respect to the incident direction of the light from 
the fiber (approximately normal to the edge of the wafer). The ridge waveguide was made 2.5 mm 
long, so that the random waveguide structure is far from the direct path of the stray light. In 
addition, uniform illumination of the front surface of the random structure inside the waveguide 
was ensured by positioning the tapered fiber approximately at the center of the input facet of the 
ridge waveguide. The spatial distribution of light intensity over the sample was imaged by an 
objective lens onto an IR CCD camera [not shown in Fig. SI (a)]. 
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Figure S 1 : Optical measurement setup, a, Schematic of experimental setup for measuring light 
transport inside the random waveguide, b, Photograph of the experimental setup, c, Schematic 
of the sample layout showing the ridge waveguides coupling the probe light from the edge of 
the wafer to the random waveguides with photonic crystal sidewalls. d, Layout of the fabricated 
structures studied experimentally. 
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